Effects of the ^-factor in semi-classical kinetic plasma theory 
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A kinetic theory for spin plasmas is put forward, generalizing those of previous authors. In the 
model, the ordinary phase space is extended to include the spin degrees of freedom. Together with 
Maxwell's equations, the system is shown to be energy conserving. Analysing the linear properties, 
it is found that new types of wave-particle resonances are possible, that depend directly on the 
anomalous magnetic moment of the electron. As a result new wave modes, not present in the 
absence of spin, appear. The implications of our results are discussed. 
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In recent years there has been a rapidly growing in- 
terest in the quantum properties of plasmas 0, 0, [J 0, 
0> 0, 3 • This has been motivated by applications in, 
for example, plasmonics quantum wells [n[ and 

ultracold plasmas ■ Common to such applications are 
rather extreme parameters, compared to most laboratory 
and space plasmas. In particular the plasma densities are 
considered to be very high and/or the temperatures are 
correspondingly low. For astrophysical plasmas, it is also 
known that strong magnetic fields (l3j may lead to vari- 
ous quantum effects being important. However, a recent 
work [8( show that the spin properties of electrons can be 
important even in high temperature plasmas of modest 
density and magnetic field strength. 

In the present Letter we will put forward a more elab- 
orate kinetic model, where the electrons are described 
using a distribution function in an extended phase space, 
including also variables due to the spin orientation. This 
model is an extension of a kinetic model used by Refs. 
14, lj|, where we here also include the magnetic dipole 
force associated with the spin. As a consequence, the 
magnetic dipole energy also contributes to the energy 
conservation law which is derived from Maxwell's equa- 
tions combined with the spin-kinetic model. This sys- 
tem is then used to study linear waves in a homogenous 
magnetized plasma. The analysis shows that the inclu- 
sion of spin gives rise to new phenomena not present 
within the usual Vlasov model. This also holds for the 
high temperature regime, where quantum effects are nor- 
mally suppressed. It turns out that effects of the anoma- 
lous magnetic moment are crucial. In particular, new 
wave modes appear with frequencies ui ~ (#/2 — V)w c , 
where g ~ 2.002319 is the electron spin gr-factor, and oj c 
is the electron cyclotron frequency. Furthermore, new 
types of wave particle interaction can take place that in- 
volve the electron spin state. We stress that none of 
these effects can be seen within quantum fluid models 

Here we are interested in effects due to the electron 
spin that may survive even when the macroscopic varia- 
tions occur on a scale longer than the thermal de Broglie 
wavelength, which is a scale normally expected to imply 



classical behavior [161 ]. In particular, we will by spin here 
mean the semi-classical properties of the electron due to 
its magnetic moment, and thus will not take into account, 
e.g., commutation relations. As a starting point, let us 
consider the evolution equations for momentum and spin 
resulting from the Pauli Hamiltonian 

H = -(h 2 /2m 2 ) [V - {iqA/h)f + ^ e cr • B + g$, (1) 

where q — — e = — |e| and m are the electron charge 
and mass respectively, A and $ are the vector and scalar 
potential, 2irh is the Planck constant, /i e = — (<?/4)eft/m 
is the electron magnetic moment, and a is the vector 
consisting of the Pauli spin matrices. Using dF/dt = 
dF/dt+(l/ih) [F, H), where F is any operator, we obtain 
in the Heisenberg picture 



V = (q/m) (E + v x B) 
S = (2/i e /?l) (s X B) , 



(2(i e /mh) V (s • B) , (2a) 
(2b) 



where s is the spin operator, E and B are the electric and 
magnetic fields, v = x, and the over dot denotes the total 
time derivative. Next, letting the number of particles 
with the expectancy value of the velocity between v and 
v+dv and spin vector between s and s+ds be given by 
dN = fdvds, we search for an evolution equation for 
/. Noting that V s • (s x B) = 0, particle conservation 
implies that the phase space density is conserved along 
fluid elements propagating in the extended phase space, 
i.e., /(r,v,s,t) = dtf + i- V/ + v • V„/ + s • V s / = 0. 
Thus, with Eqs. (2), 



dtf + v-Vf- 

2fj, e 



^(E + vxB)A (s-B) 

m mn 



~p(sxB)-V,/ = 0. 
n 



VJ 
(3) 



For the concept of a distribution function in phase space 
to be well defined, the derealization of the wave func- 
tion of the particles cannot be too large. To establish 
the limit of validity, we therefore would like to connect 
<j3j> directly to the Pauli equation. Studying the one- 
particle equation ihdt^a = H^a, where ^f a is the elec- 
tron wave function and the subscript a is a particle index, 



2 



the Pauli equation can be transformed into fluid like vari- 
ables by the transformation ty a = ^/nZ exp(iS a /fi)(p a , 
where ip a is a unit normalized two-spinor. Defining 
v a = (f /m)(yS a — ihip^Vifa) — (q/m)A and s a — 
(H/2)(pl c (Tip a , the Pauli equation can be rewritten as evo- 
lution equations for n a , v a and s Q 5], which have a 
fluid like form, although we are still dealing with sin- 
gle particle equations. These equations resemble Eqs. 
(|2"aj) and (|2b[) above, but contain several other terms due 
to the difference of a two-spinor from a classical par- 
ticle with a magnetic moment. In particular, a force- 
like term that arises in this way, even in the absence of 
spin, is the so called Bohm potential which results in 
a force Fb = (H 2 /2m)V[(V 2 y/n^) / ^/n^]. Assuming a 
characteristic (thermal) energy of the particles, the ef- 
fect of the Bohm potential is seen to be small provided 
the gradient scale length is longer than the thermal de 
Broglie wavelength [16J. Similarly, quantum corrections 
to the Vlasov equation due to wave function dispersion 
[ill , has been shown to be small whenever spatial gradi- 
ents are long compared to this quantum scale. Further- 
more, the spin properties induces higher order force terms 
F 3 = {l/m)d l [(djS k )(diS k )] [13 that are small (com- 
pared to the convective derivative), provided the gradient 
scale length is longer than the thermal de Broglie wave- 
length. Finally, higher order terms in the spin evolution 
equation (l7j resulting from the Pauli equation, can be 
neglected when the same condition is fulfilled. The ki- 
netic equation put forward here agrees with the one used 



byRefs. [Mlli, 

except that their equation did not con- 



tain the magnetic dipole force. This term is often small 
as compared to the Lorentz force, but we will show that 
it can still be important. 

A complete model is formed by combining Eq. ([3]) with 
Maxwells equations, where the current density is 



j = jfrco + VxM = 



vfidfl + x / s./>/<2 



(4) 

The sum is over particle species i with charge and 
the last term is the magnetization current due to the 
spin. Normally the spin contribution from the ions can 
be neglected compared to that of the electrons, due to 
their smaller magnetic moment. In what follows, we will 
therefore only include the electron physics, and drop the 
sum over species. The integration d£l is made over three 
velocity variables and two spin degrees of freedom. Since 
the spin vector is constructed from the expectancy value 
of the spin operator, it has a fixed length |s| = H/2. The 
spin orientation will be described using spherical coordi- 
nates. From this model it is straightforward to show that 
the energy conservation law dtW + V • P = is fulfilled, 
with the energy density W and energy flux P given by 
W = (e E 2 + ^B 2 ) /2-B-M+ / {mv 2 /2)fdn and P = 
Ex x B - M)+J [(mv 2 /2)v - (2/i e /ft)v(B ■ s))] fdVt, 
respectively. The last term in the energy flux expression 



represents the convection of magnetic dipole energy. 

Next we will use Eq. ([3]) to study linear waves in 
a magnetized plasma. Dividing the variables as / = 
/o(v, s) + /x(r, t,v, s), and B = B + Bi(r,i), the lin- 
earized Vlasov equation can be written 



d t +vV + ^(vx B ) • V„ + % (s x B ) • V, 

m n 



--(E + vxBi).V„/ 



2^e 

h 



V(s-Bi) 



V„ + (sx Bi) • V s 



fo- 



il = 



(5) 



Before proceeding further, we will for definiteness specify 
the unperturbed equilibrium distribution /q. In thermo- 
dynamic equilibrium and for a large chemical potential 
(which applies for n\ dB ^> 1, where AdB is the thermal 
de Broglie wavelength, and n is the electron number den- 
sity), the Fermi-Dirac equilibrium distribution reduces to 



no / m 
Air \2Trk B T 



3/2 

exp 



(mv 2 /2 + 2^ e s ■ B /H) 



k R T 



(6) 

This results in a zero order magnetization Mo = 
nofJ-eV i^eBo/kBT), where rj is the Langevin function, T 
is the temperature, hq the equilibrium density and ks 
is the Boltzmann constant. Next, letting Bo = -Bqz, in- 
troducing cylindrical coordinates (v±,ip v ,v z ) in velocity 
space and spherical coordinates (tp s ,9 s ) in spin space, 
noting that the term (v x Bi) ■ V„/o m © can be 
dropped when the unperturbed distribution function is 
Maxwellian, and Fourier analyzing, Eq. ([5]) is written 



i (u> — k • v) +lu ( 



d 

dtp v 



d 

' dip s 



fx 



^E + ^V(s.B i: 

m ma 



2fi c 



(sxBi)-V fl / (7) 



where we have introduced the frequencies lu c = qB^/m 
and Lo cg = 2/i e Bo/h. Note that lu c < and that uj cg = 
(g/2)uj c . Eq. ([7]) can be solved by an expansion of /i in 
the eigenfunctions ip a ((p v ,Vx) = (27r) -1 / 2 exp[— i(aip v — 
k±vx sm.(p v /ui c )]. Thus, we let 

fi = *^2gab{v±,v z ,d s )4> a (tp Vl v 1 _)exp(-ibip s ) (8) 



where a = 0, ±1,±2, ... and b ■ 
thogonality properties J 27r ^ a ^dip v 



1,0,1. Using the or- 
8 a b we find 



i (lo 



auj c 



bui cg ) g ab = Iab(v±,v z ,8 s ) (9) 



with 



lab =, 



27r r 2n 



i.E+^V(8.B 1 ) 

m ma 



V„/ 



2[i e 



H — — (sxBi) • V s / } ip* exp(ibip s )dip v dtp s . (10) 
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A useful relation when trying to write results in a more 
explicit form is the Bessel-expansion 



ipa{<Pv,v±) = ~K=y2Jb 
V2ir b 



k±v\ 



exp[i(6 - a)tp v }. (11) 



Here it is seen that the results are much simplified in the 
limit where k±vth/^c is small (where we can estimate v± 
with the thermal velocity Vth), but as is well-known (see 
e.g. Ref. [18(1), in general the conductivity tensor compo- 
nents turns into sums over Bessel functions . The conduc- 
tivity tensor o*i , defined by j l — o % i Ej , is found from Eq. 
© by expressing the magnetic field in terms of E and 
then solving for /j in terms of E using the eigenfunctions 
as outlined above. It is illustrative to divide the conduc- 
tivity tensor in two contributions a lJ — oil + crll arr „ 

^ ncc mdgii 

due to the free current and the magnetization current. 
Furthermore o~l could be divided further into contribu- 
tions due to the Lorentz force and contributions in I a b 
that contains spin which give the classical and spin parts 
of erf . However, we will not present the general expres- 
sions for cr y here, as these results are complicated and 
needs extensive analysis for a useful interpretation. In- 
stead we focus on special cases in order to point out two 
of the main new features resulting from the spin. 

Firstly, the factor (ui — k z v z — auj c — bu> cg ) in © re- 
veals that the standard wave-particle resonances are ex- 
tended to involve the spin- forces. The spectrum of res- 
onances is thereby much extended. In particular the 
combination a — 1,6 = —1 give a resonant veloc- 
ity at v z = (uj — Auj c )/k z , where we have introduced 
Aw c = (g/2 — l)eBo/w. The physics of this resonance 
is slightly more complicated than the well-known reso- 
nances with 6 = 0. In particular, it does not occur for 
strictly parallel propagation to Bo, since the coefficient 
Ii-i in (|10| becomes zero in that limit. Furthermore, 
inclusion of the spatial variations of the wave perpendic- 
ular to the magnetic field must include a finite argument 
of the Bessel functions in (fTT|) . This means that the spin- 
resonances depend on finite Larmor radius effects. More- 
over, integration over (p s kills any contribution from this 
resonance in the free current density, and thus only the 
magnetization current survives. Still it is clear that such 
resonances can be the dominating wave-particle damping 
mechanism, in case the wave frequency is of the order of 
a few thousands of the cyclotron frequency. 

The second point to be pursued in more detail is that 
the spin gives rise to new wave modes not present in 
the standard Vlasov picture. Firstly, we assume that 
the ions are immobile, constituting a neutralizing back- 
ground. Next, we consider the limit of exactly perpen- 
dicular propagation, and thus let k = fcj_x^ In this case 
modes with the approximate polarization E = E z z, and 
B = B y y = (k±E z /uj)y are possible, provided a xy and 
a xz are sufficiently small, which can be verified a pos- 
teriori (see, e.g., Ref. 18]). From Eqs. ©, © and 



(110]) the distribution function is expressed in terms of 
the electric field, that together with © determines the 
conductivity tensor components. For the given polariza- 
tion, only the a zz component is needed. Furthermore, 
for uj <C |cj c |, terms with coefficients a ^ 6 in © will be 
smaller than those with a = b, and thus only the latter 
are kept. Carrying out the (p s and ip v integrals [using 
(|11|)]. the ^-component of Ampere's law gives 



2 ; 2 2 i 2 

uj = k c + uj r 



Jo (k±v±/uj c ) 



k 2 h 2 Auj c sin 2 6 S 
4m(u — Auj c )kBT 



[J? (fcjUj_/w c )] }fodn, (12) 



where uj p = (n^e 2 /earn) 1 / 2 is the plasma frequency. If we 
drop the term proportional to J 2 due to the spin in (|12p , 
we have the usual ordinary mode, where higher terms in 
the sum over Bessel functions have been dropped, due to 
the condition uj <C |w c |. Typically the second term due 
to the spin is a small correction when the wavelength is 
longer than the thermal de Broglie wavelength. How- 
ever, it can clearly be dominant for frequencies close to 
the resonance uj « Aw c . In Fig 1 numerical solutions of 
the dispersion relation is plotted for different parameters 
£ and r\ close to this resonance, where £ = H 2 uj 2 /m 2 vf, 
T] = uj 2 c 2 /uj 2 v 2 and vt = iksT/m) 1 ^ 2 . To some extent 
these new solutions resembles the well-known Bernstein 
modes (see e.g. Ref. [Uj]). However, unless the ratio 
is larger than unity (corresponding to a high density 
low temperature plasma), the deviation from the reso- 
nant frequency will be small. For an unperturbed distri- 
bution different from the thermodynamic equilibrium ex- 
pression the shift from the precise resonance can be much 
enhanced, since the two last terms in (jTTJJ) contributing 
to almost cancel (the terms match as g/2 — 1) for 
/o given by ([6]) as considered here, but not for a gen- 
eral equilibrium distribution. Furthermore, it should be 
noted that this approximate cancellation in (|10p does not 
occur for general angles of propagation. For arbitrary an- 
gles of propagation, the denominator uj — Auj c of our spin 
term is replaced by uj — k z v z — Auj c . It is interesting to 
note that the sign of this contribution is determined by 
the last term in (|10p . In particular, in case the unper- 
turbed distribution / has a larger fraction of particles 
in the higher energy spin state, dfo/ds z changes sign, 
which should consequently change the sign of the imagi- 
nary part of the dispersion relation, leading to instability 
rather than wave-particle-spin damping. However, more 
analysis is needed to definitely confirm this conjecture. 

In the present Letter we have put forward a Vlasov- 
type equation, in a phase space extended to include the 
spin degrees of freedom. Thi s eq uation extends a previ- 
ous equation due to Refs 3, Hf, by including the mag- 
netic dipole force associated with the spin, and can be 
viewed as a semi-classical limit of the Pauli equation. In- 
cluding the spin- magnetization current in Maxwells equa- 
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FIG. 1: The normalized frequency Q = lo/Auj c plotted 
around Au c as a function of the normalized wavenumber 
K = kvt/uj c for different parameter values of the parameters 
£ = h 2 u! 2 /m 2 vf and r/ = ui 2 c 2 /uJpV 2 . The upper panel rep- 
resents possible laboratory values (for example, So = 6.0T, 
T = 10 4 K, and n = 10 22 cm" 3 correspond to ( = 1CT 6 and 
r/ = 10~ 2 , relevant for laser-plasma interaction experiments), 
while the lower are representative for extreme astrophysical 
environments (for example, Bo = 6.0 x 10 6 T, T = 10 10 K, and 
n — 3 x 10 21 cm -3 correspond to f = 1 and Tj = 10~ 3 , relevant 
for a thick hot accretion disk surrounding a pulsar). 



tions, the resulting system is shown to be energy conserv- 
ing. To gain some understanding of the model, we have 
outlined the theory for linear wave propagation in mag- 
netized plasmas and demonstrated the appearance of new 
wave modes due to the spin. These wave modes depend 
on resonances associated with both the orbital and spin 
gyration, and is much different from the well-known spin 
waves in ferromagnetic materials (l9| . 

It is of interest to discuss the connection between our 
spin-kinetic model and the ordinary Vlasov equation. Re- 
lating E and B through Faraday's law, we see that the 
relative strength of the spin force is k 2 h/muj. This pa- 
rameter is often small, which explains why the classical 
Vlasov equation in many cases is a very good approxi- 
mation. On the other hand, using this parameter it is 
easy to underestimate the significance of the spin force. 
Firstly, wave-particle resonances and/or new wave modes 
associated with resonances, as the one described above, 
may enhance the significance of spin effects. Secondly, 
when the electric field is perpendicular to Bo and the 
wave frequencies are low, the E x B-motion of electrons 
and ions give a zero current to leading order in uj/oj c , 
whereas similar cancellations does not occur for the mag- 
netic dipole force. Furthermore, even if the electric field 
is perpendicular to Bo, the dipole force may have a paral- 
lel component leading to a large current in this direction. 

The exploration of Eq. {3]) is a very rich problem. Be- 
sides the more obvious generalization to a complete linear 
theory, including arbitrary angles of propagation, there is 
a large set of nonlinear problems that should be studied. 
The results by Ref. [8| indicate that the significance of 
the spin force can be increased in the nonlinear regime as 



compared to the linear one, when the spin dependence of 
the ponderomotive force can separate particles of differ- 
ent spin. Finally, it should be stressed that the current 
model is intended for the " weak quantum" regime where 
the characteristic length scale is longer than the thermal 
de Broglie wavelength, and the Zeeman energy is smaller 
than the thermal energy. A very interesting problem is 
generalizations to the regime where the Zeeman energy is 
comparable or larger than the thermal energy. For such 
parameters, found in astrophysical settings, effects such 
as Landau quantization [20j will enter the picture. 
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